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ABSTPAf  T 


The  mathemat ica 1  modeling  of  combat  engagements  by  heterogeneous  tactical 

forces  is  discussed  and  a  new  superiority  parameter  is  introduced  which 

represents  a  measure  of  the  effectiveness  of  the  opposing  forces.  This 

parameter  is  a  function  of  both  the  Quantitative  and  qualitative  strengths  of 

the  opposing  forces.  It  is  defined  as  a  product  of  the  transpose  of  the  left 

dominant  eigenvector  of  the  governing  matrix  of  attrition  coefficients  and  the 

force  strength  vector.  This  parameter  is  used  to  establish  superiority 

ci iteria  for  combat  engagements.  The  application  of  this  new  concept  to 

heterogeneous  forces  is  illustrated  for  the  case  of  a  "one-nn-twori  tactical 

i  _  r 

engagement.  r 


1  ■  I rtroduc.  t  i  on 


the  Lanchester  equations  fo  the  lie* et ogeneous  forces  engaged  in  a  direct 
fire  to mb at, can  he  represented  symbol irall,  by 


whe  r  e 


x  --  <;x 


X  ;M  M 


ard  IS  the  nati  i>.  nf  attrition 


C  DPf  ♦  1  ■'  1  ei '  t  a  . 


He.  e  M 


N 


;  o  1  umi 


matrices  denoting  the  force  (weapon)  numbers  of  the  Blue  and  Red  forces 
respectively.  The  braces  are  used  for  convtiiienrn  throughout  th:s  report  to 
denote  column  matrices.  As  discussed  by  Taylor*.  several  authors  have 
introduced  aggregated  force  parameters  to  obtain  a  better  insight  into  the 
effects  of  varying  force  size,  the  weaoon  allocation,  and  the  we  at' or  effective- 
firing  rate  represented  by  Eq .  (1/.  This  aDbroarh  im-m.-pri  essential), 
computation  of  two  quantities:  the  aqareaated  force  strengths  F  and  F  fr- 

-  -  B  R 

the  Blue  and  Red  forces,  respectively,  which  may  be  defined  as 

~  =  v1  H  3 

B 

F  =  w  rN  ( -•  1 

R 

where  v  and  w  are  positive  vectors  representing  relative  .aluc-s  of  individual 
weapons  while  H  and  N  are  column  matrices  representing  the  force  (weapon) 
numbers  for  the  Blue  and  Red  forces,  respectively  Starting  with  the 
assumptions  (3)  and  (A)  ,  a  number  of  authors  used  the  eigenvector  method  of 
computing  v  and  w  (see  Section  3).  In  this  method,  however,  both  vectors 
could  only  be  determined  up  to  a  constant  multiplier.  A  simple  sr a li nq 
relationship  was  proposed  b  /  Dare  and  James  as  £  v  =  £  w  --  1.  Other 

plausible  scaling  methods  were  introduced  by  different  authors.  For  example. 
Taylor1  in  his  boot  on  Lanchester  models  of  warfare,  mentioned  three  othr-' 
sc  hemes’*  °  for  cdmputinq  the  scaling  factors  for  the  weapon  .alue  vecto  s  v 
and  w 

This  report  introduces  a  new  wav  of  detei  mining  the  ei  cier.vec  to»  s  v  -.n.J  w 
through  a  mathematical  represent  at.  i  on  of  the  complete  anal/tical  sc  lot  ion  and 
the  application  of  the  dominant  lef1-  eigervec  tor  o'  the  matri-  ;!ac 


T 

approach  leads  also  to  a  new  measure  of  effectiveness  u  X  or  its 

to 

corresponding  nond 1  mens i ona 1  superiority  parameter  S  ,  both  of  which  determine 

a  priori  the  direction  of  the  outcome  of  the  battle  engagement  without  the 

necessity  of  solving  the  differential  equations  representing  the  engagement. 

The  application  of  this  new  method  cf  analysis  is  demonstrated  here  for  the 

case  of  a  tactical  engagement  involving  three  different  weapons:  one  Elue  and 

two  Red  weapons.  In  order  to  illustrate  the  nature  of  the  analytical  solution 

and  ,  in  particular,  the  presence  of  a  dominant  term  in  the  solution  which  is 

used  to  derive  the  superiority  parameter  5  ,  a  complete  discussion  of  the 

o 

general  analytical  solution  is  included  in  the  Appendix. 


r .  Lanchester  Equa t i ons  for  Heterogeneous  Force  Combat  :  "  Many-on-Many  " 

Combat  Engagements  . 

The  two  separate  governing  equations  for  the  heterogeneous  force  combat 
engagements,  from  Eq  .  <  1 )  »  are  given  by 


dM/dt  = 

-  BN 

( 5 

dN/dt  = 

-  AM 

( h 

where  M  and  N  are  the  matrices  representing  the  Clue  and  Ped  forces  icn 
weapons),  respectively.  The  mathematical  model  for  such  engagements  is 
represented  as  the  state  variables  M  and  N  interacting  with  each  other  through 
the  appropriate  attrition  coefficients  A  and  B  of  A  and  B  as  shown  in  Fig. 

■  I  !-  i 

1  Here  the  attrition  coefficients  tat-e  into  account  the  f  i  r  i  rig  rate 
effectiveness  and  the  allocation  o-r  weapons.  The  matrices  A  and  B  are 
nonneqat i ve ,  a  propet  tv  which  will  be  used  in  the  subsequent  aaalvsis.  Ah  r 


3 


M  ,  N  >  0.  f Cm  al  1  i  >■'> 

and  when  this  condition  15  not  sat is'isd)  i.e.  when  c^e  of  the  forces  is 
annihilated  ( when  a  particular  m  or  M  =  0),  Eqs.  (5<  arid  ( 5)  must  be 
reconfigured  by  striking  out  the  appropriate  rows  and  columns. 


Eguations  (5)  and  16)  can  be  conbined  mtu  a  ample  matrix  equatic; 

dX/dt  =  X  =  CX  f  1  '• 

where 

X  =  ( M  N  •  2) 

and 


O  -B 
-A  O 


(Si 


If  the  i  t.h  row  in  f.  contains  all  zeros  this  means  that  the  corresponding  X  =  0 
and  the  X  -force  is  not  being  fired  at.  If  the  jth  column  in  C.  contains  all 

zeros  then  the  X -force  is  being  inactive. 

1 

The  left  handed  eigenvector  q  of  the  matrix  C  arid  its  nor  respond  1  no 

eigenvalue  are  determined  from 

CTq  -  >  q 


The  eigenvector  q  can  be  pa1  titioned  into  rows  corresponding  to  the  Blue  and 
Fed  forces,  respectively,  so  that 


H 


ft 
’  1 


(C  V 


1  w 

?. 


where  v  and  w  are  the  column  vectorb  wmie  c  and  i:  are  constant  multipliers. 

1  2 


It  should  be  observed  that  Eq .  (1)  must  represent  a  set  of  oupled 
equations.  Uncoupled  sets  are  not  admissible  in  the  present  analysis  because 
they  represent  unrelated  (uncoupled'  tactical  engagements .  Examples  of 
uncoupled  and  coupled  engagements  are  shown  in  Fig.  2.  where  arrows  indicate 
the  direction  of  fire.  For  these  examples.  the  corresponding  matrices  of 
attrition  coefficients  are  shown  by  Eqs.  (11a)  and  (lib)  where  x’s  indicate 
the  nonzero  coefficients  in  C. 


-  v.  u 
0  -x 
0  0 
0  0 


0  -x 

0  0 


.uncoupled  sets 


(Ha) 


.coupled  sets 


(lib) 


The  test  whether  the  matrix  C.  represents  uncoupled  sets  of  equations  is  simply 
whether  it  can  be  reduced  tG  two  or  more  nonzero  diagonal  submatrices  bv 
rearranging  rows  and  columns  in  C,  while  all  other  off-diagonal  submatrices 
are  zero.  Otherwise,  the  matrix  C  represents  coupled  sets  of  equations. 


Substituting  Eqs.  (8>  and  <!0)  into  Eq.  <9), 


O  -A  c  v 

l 

-BT  O  c  w 


from  which 


A  c  w  =  ‘  r  v 
2  1 


BT  v  = 

i 


■\  r  w 

2 


T 

Premultiplying  Eq.<13)  by  B  and  then  substituting  Eq.ila)  and  pi  emu  1 1 1 p 1  ,  i  ng 

T 

Eq.(lA)  by  A  and  then  substituting  Eq.(i3),  the  following  pair  of  equations 
is  obtained 


T  T 

(A  B  - 

>  2 1  >c  v  =  0 

1 

t  15 

T  T 

BA  - 

\2I  )  c  w  =  0 

(  16 

2 


If  the  dimension 

s  u(  A  are  nxm  and  those  of 

B  are 

my.  n 

and  1 

f  m>n,  then 

there 

will 

be  n  common 

2 

eigenvalues  K  =  u  in  Eqs. 

(  15) 

and  ( 

16). 

I n  add  1 1 1  on . 

there 

wi  1 1 

be  (m~n> 

eigenvalues  of  Eq .  (15) 

equa  1 

to 

zero 

for  which 

thei  r 

corresponding  eigenvectors  can  be  obtained  from 

T  T 

A  B  c  v  =  O  (17i 

i 

If  those  eigenvalues  were  not  equal  to  zero,  this  would  then  imply  that  in 
addition  to  the  2n  eigenvalues  \  =  ±  Vij  in  C  there  would  be  additional  2(..r-n! 

eigenvalues  for  a  total  2n+c?(m-n>  =  2m  which  would  be  greater  than  the 

required  number  <m+n). 

T  T  T  T 

The  matrix,  products  A  B  and  B  A  in  Eqs.  (15)  and  (16)  are  nonneqative 

since  all  A  and  B  are  either  positive  or  zero.  Therefore,  according  to  the 

M  L  i 

Frobenius-Perron  theorem  <'Ref .  17 ,  p.  193T,  Fqs.  (15)  and  (16)  have  a  real 
nonneqative  eigenvalue  \  which  is  either  equal  to  or  exceeds  the  moduli  of 
an\  other  eigenvalue.  Also  to  this  maximal  eigenvalue  there  co  responds  an 
eigenvector  with  all  nonnegative  elements  which  will  be  denoted  as  c^v  and 

c  vy  for  Eqs.  <  15)  and  (16),  respectively.  Since  c  v  and  c  are 

21  1121 

submatrires  of  the  left  picienvector  of  C  it  follows  from  Fa.  1  10)  -hat 

u  =  (  c  v  c  w  }  >  1 9  • 

1  112  1 


6 


which  can  be  determined  directly  from  cq.  (9).  When  commuting  v  and  w  from 
Eqs.  (15)  and  (lhi  any  suitable  scaling  method  can  be  used.  For  example  the 
largest  element  can  be  made  equal  to  unity*  but  the  actual  scaling  between  v 

1 

and  w  is  acco-T.pl  i  shed  throuqh  the  factors  c  and  c  to  be  determined  latei  . 

1  12 

Transposing  Eqs.  (13)  and  (141  with  v  =  v  and  w  =  w 

i  1 


c 


T 

W 


A 


2  1 


A  c  v 

i  i  i 


i  j9> 


T  T 

-  c  v  B  =  \  c  w 

11  12  1 


t  20 ' 


Next*  premultiplying  the  above  equations  by 


e  =  {  1  1  ...  1  } 


(El  ) 


and  eliminating  the  following  relationship  between  r  and  c  is  obtained: 


f  c 

E 

T  T 

'  v  Do  v  o  N 

2 

i  i 

— 

— _ _ _  _ 

r 

T  4  T 

.  J 

w  Ao  w  e 

-  i  i 

i  i 


c 

v  Be  v  e 

2 

i  i 

—  - 

— -  - 

C 

T  .  T 

1 

w  Ae  w  e 

i  «  7 

i  i 


The  neqative  sign  in  Eq.  (23)  is  selected  to  satisfy  Eqs.  (19)  and  (20'  in 

which  A,  R,  v  ,  and  w  are  all  nonneqative  matrices  and  *  is  a  positive 

ii  i 

number.  Co  incident lv*  Eg.  (23)  is  of  the  same  form  sugoested  b.  Ta.lor 
without  proof  for  the  so-called  "summed  result"  in  interpreting  the  aggregate 
force  strengths.  Equation  (23)  provides  a  rational  relationship  between  v 

i 


7 


and  w  when  these  eigenvectors  are  calculated  fro"-  Cos.  .  !c  '  a  d 

i 

It  should  be  noted  that  i  r-  practice  the  left  dominant-  eiae  .  a  1  ~e  •  a 

i 

its  corresponding  eigenvector  q  can  he  obtain,  d  direcM-,  f  -  r  :  t,  ~q.  ■  at'  er 
than  from  Eqs.  !15>,  il6),  and  •  ?3  >  .  Computer  programs  such  a  s  MATLAB  are 
available  which  can  be  conveniently  used  to  find  the  dominant  e  ieen  .  a  I  v.e 

i 

and  the  dominant  left  eigenvector  q  . 

Because  of  the  relationship  (A.  <?8).  the  general  l  red  left  eiaen.e.  t.;rc  q 
are  orthogonal  to  the  right  eigenvectors  p.  Tins  orthogonalit  c roper t» 
expressed  as 


r 

q  p  =  1  fo>  l  i 

>■  .1 

-  0  for  i  -r  j 


J 


M  ’ 


This  property  is  used  next  to  derive  a  scalar  parameter  from  the  relation  of 
the  governing  equations  by  premu 1 1 ip  1  /  i ng  X  b,  the  transpose  of  toe  domiia,' 
left  eigenvector  q^  (Corresponding  to  the  largest  eigenvalue  ■  .  T t .e- 

resiilting  parameter  is  a  measure  of  merit  for  the  combat  engagement  . 
Consequently,  employing  Eqs.  (A.S5)  and  (A.  14b)  and  noting  that  fo>  'he 

notation  used  ir  Appendix  A  \  ='*.  1  and  q  ~q 

i  i 


T 

q  t x 


t-o 


q^  P  f  exp<  X  t  'D  J  QTX^ 


Tr  1  ,21 

q/p  p 

.  .  ]  f”  exp!  ’■  t  >.'£? 

«Tx„) 

- 

J 

O 

0 


:  l  i  > 

.  .  ]  o  o  . . .  I , 

e«p  ( 1  t  J 

o 

if  ,  T 1 

!  i  q  ix 

L 

J  1 

i 

e  v:  p  (  *■  2  t  > 

j  :  'l  :  «.) 

i  j  2  \  ! 

!  !  q  1  i 

:  !  1  ! 

o 


T  T 


e  -  c  : 


[  1  0  0  ...  3  O  O 


I 

1  2 

.  .  .  q  q 


X  ! 

<->  I 


lqT  X  1  *  qTX 

L  1  °J.-„ 


Sut-st  i  tut  nq  now  Ea  .  18.  tor  q  into  Eq.  ( £?5  >  and  observiuq  that  X  =  <  M  N  ) 

1  •  0  1)0 


it  follows  that 


1' 

q  X 
l  o 


r  T  T  1 

lev  c  v*  1 

11  2  1 


M 

< 

N 


c  j  v  rM  ‘  •:  -  c  ■  wTN  1 

1  [  1  o  2  l  1 


rD  T 

v  Bo  v  t*  ..  i  2 


=  c 


pH  r  —  1  l  «tn  ] 

»  {  »  o  l  T  T  )  1  O 

'■  w  Ae»  v*  c>  ^ 


wttrp  for  convenience  c  can  be  taker,  as  one.  Ber  ause  v  and  w 

t  i  i 


both 


T  T 

ncroeQat i ve.  the  terme  v  M  and  w  N  a. e  both  positive.  The  scalar  quant i t , 

1  O  1  O 

T 

q  X  derived  from  the  predominant  term  in  the  solution  represents  the 
difference  between  aagregate  Blue  and  Bed  feces  and  can  therefore  be  ta*-en  as 
a  measure  of  merit,  i.e.,  a  superiority  indicate'  for  the  opposing  forces.  It 
should  be  noted  from  the  above  that  the-  anal',  t  ical  solution  for  X  is  a  sc  ic  1 

t  /  wh.  ,  ,  r  i  a  t'.e  anqeo  i  r- 


{if  t.e-  ms  with  mu  1  1 1  p>  1  v  t  nq  factors  e-p-’  ! 

V  '  1  2 

descending  or'der  for  *  as  ‘  ,  .  .  .  f> .  >  ’ .  o .  .  .  .  -  * 

term  with  *•  •  1  o  t  tie  dominant  term  in  the  solution.  1  1 


Thus  the 
o  a  repea -ed 
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eigenvalue*  or  if  other  eigenvalues  are  close  to  the  dominant  eigenvalue. 


use  of  q  X  as  a  measure  of  merit  may  not  be  appropriate.  Consequent  i  »  •  |f  >  ~- 

1  o 

important  to  check  the  distribution  of  eigenvalues  before  aop i . : on  t*-i ; 

measure  of  merit.  In  most  practical  cases,  however.  the  '  1  ’s  are  we  1 

separated.  For  example,  for  equal  coefficients  in  both  A  and  B  ;i.e.  A  =- 

and  B  -b  for  all  i  and  j  where  a  and  b  are  constants'  the  eicen\ aloes  ac  i 
1 1 


*  *0*0,... *0*0, 


=Vnvr =b  and 


■  /iiD  i  ,d 


-Vmb/na.  .  .n-tunwl  where  m  and  n  refer  to  the  number  of  Bl  ue  and  Red  forces . 
respectively  (see  Appendix).  It  is  interesting  to  obsei ve  that  fci  r> e 
multiple  eigenvalue  cases,  the  dominant  eigenvector  q  contains  :ero  elements 

T 

which  means  that  some  force  strengths  are  not  included  in  q^X  . 

To  reduce  Eq .  (£6)  into  a  nondimensional  parameter,  the  eouation  can  be 

T 

normal  iced  with  respect  to  c  v  M  by  introducing  a  nondimensional  superior]!, 

1  t  o 

parameter  5  such  that 
o 


5  =  qTX  c  VTM 
i  »  l  o 


where 


qTX 
1  o 

c  vT  M 
1  1  o 


v  Be 

f  _i _ 

T 

v  e 
i 

T 

V  2  W’N 

I  1  A 

^  w  Ae 
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T 

W  <3 

1 

J  VTM 

1  1 

-  1  - 

i-'1 

o 

Ta 

w  Ae 
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T 

i 

l  2  V1  M 

vTBe 
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T 

V 

1 

1  T 

w  N 

1  13 

<■;  c. 

Di  -J  1 

anti  T’ 

at  e  a  1  1  1 

The  nondimensional  parameters  S ,  S  ,  and  T  are  all  1  ndenenclent  of  the  ihont 


of  The  arbitrary  multiplvinq  factors  witt  and 


The  i: up i 1  ■  i  c  i‘ 


parameter  5  is  a  scalar  quantity  measuring  the  contributions  of  tne  dominant 


term  in  the  solution  for  the  force  levels  at  any  given  time.  It  should  also 
be  noted  that 


0 

:  5  < 

o 

1 

f  or 

i  •  ■& 

o 

5  = 

o 

0 

f  or 

ii 

0 

-cv> 

•  5 

o 

o 

f  or 

0  <  $ 

a 

.  30  i 


r 

The  parameters  q  X  or  S  can  be  used  as  measures  of  effectiveness  to 
1  o  o 


determine  the  outcome  of  a  particular  tactical  engagement  without  actually 


solving  the  differential  equations  subject  to  the  restrictive  conditions 


T 

discussed  earlier  for  the  measure  of  merit  q  X  .  Thus  as  in  the  case  of  a 

1  o 


one-on-one"  engagement  the  Blues  forces  win  the  engagement  when  'I  >  1  or 


when  S  >  0*  while  the  Red  forces  win  when  $  «.  1  or  when  S  0.  When  l  -  1 

o  u  o  u 


or  S  -  0,  both  sides  reduce  their  strenqths  proportionately, 
o 


It  should  be  noted  that  as  in  the  case  of  a  "one-on-one"  engagement  the 

normalized  force  ratio  I  for  the  general  case  consists  of  a  product  of  two 

u 

components  representing  the  qualitative  and  quantitative  ratios.  The 

T  T  T  T  1  2 

Qualitative  ratio  is  C  (w  Aelfw  t?  )  /  (  v  Be  >  <  v  e/3  and  the  quant  i  t  a  t  i  v  e  ratio 

1111 

is  <  vTM  )  '  !  wrN  i . 
i  o  to 


Aqqr eg a ted  Force  51 r enqths 


to  reduce  the  mathematical  repi esenta t l on  of  the  ’  'mai  ,  r 
engagements  to  a  "one- on- one"  rppi-esentat  ion  the  aggregated  force 
may  be  introduced  into  Eqs .  (5)  and  <6  >  .  Premultiplying  F.q .  h, 
Fq.  !oi  tv  c  w *  ,  the  fnllrwinq  d i f f erent i a  1  equations  are  obtained: 

z  l 


-mar  , 


s  1 1  enq  t  hs 


r 

c  v  and 

i  i 


1 


T 

c  v  M 

i  i 


(31  > 


c  VTBN  =  P c  w '  N 
11  2  1 


c  w  N  - 

2  1 


c  w1  AM  ~  ::<c  v1  M 
2  1  11 


where  the  overdots  represent  d i f f er ent i a t i on  with  respect  to  time  and  the 
right  sides  of  these  equations,  with  positive  constants  and  /'?,  have  at  this 
point  been  introduced  arbitrarily.  Introducing  the  aggregated  force  strengths 
F  and  F  from  Eqs.  (3)  and  (4)  as 

B  R 


F  vTM  =  L  VTM 

B  11 


F  =•  wTN  =  -c  wTN 

R  2  1 

where 


v  =  c  v 

l  l 


w  =  -c  w  N 

2  i 


Hence,  substituting  Eg.  (33)  and  (34)  into  Eqs.  (31)  and  (33) 


F  -  -;>F 

B  R 


F  =  -;iF 

R  B 


which  are  the  governing  equations  for  a  "one-on-one"  engagement  between  the 
segregated  forces  F  and  F  ;  however ,  a  word  of  caution  is  appropriate  here. 

B  R 

The  solutions  for  F  and  F  are  onlv  appr o < i ma t i ons  to  the  e  act  solutions 

B  R 

because  they  do  not  incorporate  the  condition  that  one  of  the  M  or  N  can  be 

l 

reduced  to  zero  independently  of  others.  In  the  aggregated  solution  of  Eqs. 
(37)  and  f33>  all  components  of  F  or  F  are  reduced  to  zero  simultaneously. 

B  H 


Since  H  and  N  are  arbitral'/,  in  order  for  the  right  sides  of  Fqc 


and  <3S)  to  be  valid  the  following  relations  must  be  true: 


x  x 

-C  v  B  =  /3c.  w  or 

11  2  1 

T 

c  v 
i  i 

=  /?C  w 
2 

T  .  T 

-C  w  A  =  cvc  v  or 

-ATc  w 

=  etc  v 

21  11  22  11 


Solving  for  v  and  w  , 

li 


Tt  T 

A  B  c  v  =  -/?A  c  w  =  /?ac  v 

11  2  111 


(ATBT  -  af>l)c  V  =  O 

i  i 


XT  X 

B  A  c  w  =  -otB  c  v  -  cx(3c  w 

2  1  11  '  2  1 

(BTAT  -  ctpl) c  w  =  O 

2  1 

Hence,  comparing  Eqs.  (41)  and  <42)  with  Eqs.  (15)  and  (16) 

=  a/? 

Similarly,  by  comparing  Eqs.  (39)  and  (20)  and  Eqs.  (40)  and  (19) 


i 


which  is  the  same  result  derived  from  the  scaling  method  assumed  by  Holter 

.  .  ,  <5 

and  Anderson. 


The  above  analysis  has  demonstrated  that  a  unified  theory  for 
analysis  of  Lanchester  equations  for  heterogeneous  forces  engaged  in  a  di 
■tire  combat  can  be  developed.  This  is  accomplished  by  correlating 
quantities  v  w  .  ci  and  /?  ,  used  previously  by  other  authors,  to  the 


( 3Q  ) 

( 40 ) 


(41  ) 

(42) 


(  43) 


(44  ) 


the 
rec  t 
the 
left 


dominant  eigenvector  of  the  governing  matrix  of  attrition  coefficients  and  its 
dominant  eigenvalue.  The  theory  provides  for  a  proper  scaling  of  the 

T 

eigenvectors  v  and  w  and  leads  directly  to  the  measure  of  effectiveness  q  X 

’ll  i  o 

or  the  superiority  parameter  5  .  Although  the  use  of  the  left  dominant 
eigenvector  introduced  here  for  the  combat  analysis  is  new,  similar  approach 
has  been  used  in  the  past  in  other  areas.  For  example  the  use  of  the  left 
dominant  eigenvector  in  the  cohort  population  model  that  leads  to  the  total 
reproductive  value  of  the  population  as  described  by  Luenberger  CRef.7, 
p.  1853. 


A.  Numer ica 1  Example;  " One-On-Two'  Tac t i ca 1  Engagement 


The  proposed  method  of  analysis  of  heterogeneous  forces 
illustrated  here  for  the  case  of  a  tactical  engagement  of  one  Blue 
two  Red  forces  as  shown  in  Fig.  3.  For  this  example  the  equations 
the  engagement  have  been  assumed  as 


combat  is 
force  and 
descr ibing 


' 

0 

-3/2 

-1 

x  ! 
2 

-4/3 

0 

0 

X 

L  3- 

-2 

0 

0 

It  can  be  verified  easily  that  for  this  case 


X  =  (X  X 

X  }  =  (M  N  N  } 

1  2 

a  112 

q  =  (1 

-3 '4  -1/2} 

i 

T 

v  =[13; 

w'  =  n  2/31 

i 

1 

y 

11 

II 

ru 

14 


X 


and 


■$ 

o 


t 

nTTTfP'fm 

2  3 


I n  F igs . 

4,  5  and  6  the  variation  of 

X 

i 

X 

2 

and  X  with 

3 

the 

nondimensional 

time  r  is  shown 

for  X  = 
o 

{ 200 

o 

o 

o 

c 

,  (125  100  100), 

S  l  *SJ 

{ 100  100  100) 

which  correspond  to 

f  =  1.6, 

L> 

1.0, 

and  0.8, 

respect l vely . 

The 

solutions  for 

X  were  obtained  by 

numer ical 

l ntegr at  ion. 

For  Case  1  (Fig. 

4  1 

in  which  J?  = 
o 

1.6,  X  =  0  when  t 

3 

=  0.57  and  at 

that  time 

the  c ,  *■  i  o  n  of 

X 

1 

allocated  to  fight  X^  is  reallocated  against  X_^,  requiring  that  the  equations 

must  be  reformu 1  a  ted .  Subsequently  the  X^-force  is  reduced  quickly  to  zero  by 

the  X^-force.  As  the  battle  progresses,  the  superiority  parameters  S  shows  a 

rapid  increase  (see  Fig.  7  ).  Here  the  initial  positive  value  of  S  is  an 

indication  of  the  superiority  of  the  Blue  force.  For  Case  2  in  which  $  =1.0 

o 

,  the  opposing  forces  diminish  their  strengths  gradually  without  any  side 
gaining  a  clear  advantage  (see  Fig.  5).  Here  S  =  0.  Case  3  for  which  = 

0.8  represents  a  rapid  demise  of  the  X  -  force  (see  Fig.  6).  For  this  case  S 
starts  with  a  negative  value  which  diminishes  as  the  battle  progresses  (see 
Fig.  7).  Tne  initial  negative  value  of  S  is  an  indication  of  the  inferiority 
of  the  Blue  force. 


Since  the  dominant  left  eigenvector  is  the  key  parameter  the  proposed 

measure  of  merit  5  .  the  eigenvectors  q  are  shown  as  examples  in  Appendix  B 

o  i 

for  the  cases  of  ’m-on-n"  engagements  for  which  all  attrition  coefficients 
are  equa  1  . 


5 .  Conclusions 


The  premultiplication  of  the  solution  vector  X  =  iM  N)  for  the  engaging 
forces  by  the  transpose  of  the  dominant  left  eigenvector  q  of  the  matrix  of 
attrition  coefficients  helps  to  clarify  the  meaning  of  the  eigenvector  methods 


lb 


of  determining  the  relative  values  of  heterogeneous  forces  in  a  combat 

engagement.  Although,  the  eigenvector  methods  have  been  around  for  a  long 

time,  problems  existed  in  deciding  on  the  best  method  of  scaling  the  resulting 

eigenvectors.  The  present  report  demonstrates  clearly  the  meaning  of  the 

individual  eigenvectors  v  and  as  components  of  the  dominant  left 

eigenvector  of  the  governing  matrix  of  attrition  coefficients.  This  in  turn, 

allows  for  a  proper  determination  of  the  relative  scaling  of  eigenvectors  for 

the  two  opposing  forces  and  it  leads  to  the  nondimensional  superiority 

parameter  S  which  can  be  used  as  a  measure  of  merit  and  an  indicator  of  the 
o 

final  outcome  of  the  battle  engagement. 

The  practical  utility  of  the  proposed  measure  of  merit  can  be 
demonstrated  by  obtaining  numerical  solutions  to  heterogeneous  force 
engagements  while  at  the  same  time  computing  the  parameter  S  introduced  in 
this  paper.  The  sign  of  the  parameter  5  and  the  extent  to  which  it  is  either 
increasing  or  decreasing  can  be  used  as  a  measure  of  superiority  of  the  Blue 
forces  in  relation  to  the  Red  forces.  Although  the  method  is  valid  only  for 
constant  attrition  coefficients,  it  could  also  be  of  some  value  for  cases  with 
variable  coefficients  where  the  instantaneous  values  of  5  could  be  used  to 
indicate  the  general  trends  as  the  battle  engagement  progresses. 
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APPENDIX 


Analytical  So iut ion  of  Lanchester  Equat ions 

a .  General  Solution  of  X  -  CX 

The  general  solution  of  the  fundamental  equation 

X  =  CX  (AD 

can  be  obtained  as  a  product  solution  of  the  form 

X  -  Pp  (AS! 


where  P  is  the  matrix  of  general  i2ed  right  eigenvectors  and  p  is  the  column 
matrix  whose  elements  are  function  of  time.  Substituting  Eq.(A2)  into  (Al> 

X  =  Pp  =  CPp 


Hence 

where 


■f  =  P  tCPp>  =  J  p 


<  p 

:l>  <2> 

u> 

<P  ■  • 

.£*>  .  . 

.  1> 

.  <2> 

.  (V) 

% 

<P 


(A3  ) 

( A4  ) 


and 


J  =  P^CP 


r  ,<*> 


o 


o 


( Aba  ) 

( A6b  ) 


is  the  Jordan  canonical  form  of  the  matrix  C,  modified  to  take  into  the 
account  the  dimensionally  of  the  matrix  elements.  The  governing  equation 
matrix  C,  its  eigenvalues,  and  the  matrix  J  are  of  the  dimension  Dtime.  Also 
for  reasons  explained  later  J  is  taken  here  as  the  lower  triangular  matrix 
instead  of  the  conventional  upper  triangular  matrix. 


A  typical  term  in  Eg. (A3)  is  of  the  form 


•  A"  ) 
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:t  > 
Y 


It) 


(  u 


\ 

J 


f  A9  > 


where  J'1'  is  essentially  one  of  the  three  types*  or  combinations  thereof,  as 
discussed  oeiow. 8 


Type  CiP: 

JU>  --  \'1' 

when  the  algebraic  and  geometric  multiplicity  of  the  eigenvalue  K"~' 
to  one. 


i 

<  A  1 0  ) 
is  equal 


T  y  pe  C  i  i  P  : 


xu>i 


(All) 


when  the  algebraic  and  geometric  multiplicity  of  the  eigenvalue  A‘  'is  equal 
to  k  with  k>2  and  the  corresponding  eigenvectors  are  linearly  independent. 


J y pe  c  i i i J 


—  X 


0 

0 


1  0  0 
1  1  0 
0  1  1 
0  0  1 


0  0  0  0 
0  0  0  0 


0  0 
0  0 
0  0 
0  0 


0 

1 


(  A 1  2  ' 


k>:k 


when  the  algebraic  multiplicity  of  the  eigenvalue  K"  1  exceeds  its  geometric 
multiplicity.  The  expression  for  J'w  has  been  modified  here  .from  its 
standard  textbook  form  in  which  only  the  diagonal  terms  are  equal  to  '  ",  to 
ensure  proper  dimensionally  of  all  its  elements.  The  eigenvalues  *  in 
two  different  J  "’’’s  can  be  the  same.  Also  the  7'y pe  ■  i  n  ’>  mar  appear  in 
combination  with  Type-e  C  ?  and  i  i  P.  When  k=l  the  Type  i  \  ;  '  reduces  to  7  vp-  j 
C  i  P  . 
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It  is  easy  to  verify  that  the  solution  to  Eq.<A7)  with  J given  b'. 
Eq(Al?)  (i.e.  .Type  v’itO)  can  be  obtained  by  solving  first  for  p  and  then 
substituting  it  into  the  equation  for  ip"''.  This  process  is  repeated  until 
the  solution  for  1  ‘  is  found.  These  solutions  can  be  written  concisely  as 


<P 


,u  f  K'1  '  t 
=  exp  (X  t )  I  +  — ~ 


<\iL't>2 
2  1 


•  •••  K 


=  expo,  t )  1 1  + 


(- 


”■  '  H  .  h*  . 


K 


( A13a) 


i A 1 3b  ) 


.  C L  >  _<l>  ,<l> 

=  e  x  p  ( A  t )  ;D  si 


( A 1  3c  > 


where 


>!l  * t  '  \<t  >  t  )2 

iD  =  I  +  — — i  Hr'" 


1  * 


x‘l>t 


2! 


H2  r  +  o/  ‘  'O"*  Hk  l 
*  ( k - i  > »  H 


i 

u  > 


(\<0t!Z/2 
(h‘  1 '  t )  3  /  6  (X<1  t )  2/?  \tl’t  1 


<  A 1  a  ) 

(AUb  ) 


and 


'*l>  r  J  l>  -.1  l  > 

.4  =  (  .4  Sf 

t  z 


,*i0  ) 


A 1 5  ) 


is  the  column  matrix  of  constants  of  integration.  The  H-matrices  are  given  b. 


H  =  H  = 

i 


0 

1  0 

0  1  0 


'  A 1  a  > 


kxl 


H 


H2  = 


0 

0  0 
1  0  0 


o 


o 


0  0 


1- :  i 


etc  . 
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which  finally  leads  to 


Hk  =  O  <A16> 


Substituting  Eq.(A13c)  into  <A4)  and  (A2),  it  follows  that  the  general 
solution  of  X  =  CX  can  be  written  as 


X  =  Py>  =  P  f  exp  ( A  t )  ;D  J  A 

i.  A 1  9  '• 

where 

[~  exp  (At )  .73  J  =  exp  ( \!1/t )  exp  (  x'2'  t  )X>'2‘  ...  J 

i  A20  ) 

II 

* — i 

N 

< A21  ) 

At  time  t=0-  X=X  and  therefore  from  Eq.(A19> 
o 

X  =  P.tf 

(A22> 

since 

o 

exp(0)  =  l  and  ID=  I  for  t=C>.  Hence 

A  =  P_1X  =  QTX 

< A23 ) 

where 

o  o 

the  matrix  Q  is  defined  as 

-1  T 

Q  =  (P  > 

( A24  > 

Hence 

from  Eqs.(A21)  and  ( A24 ) 

X  =  P  f  exp  (At  >  ID  J  QTXo 

i  A25 ) 

If  the  matrix  C  is  not  defective-  i.e.  its  geometric  multiplicity  is  equal  to 
the  algebraic  multiplicitv  for  all  eigenvalues  >.<l  .  then  :DA,-1  and  the 
solution  for  X  simplifies  to 

X  =  P  f  exp  (At)  I  0TX  (A26> 

1  J  o 

b .  Right  and  Left  General  1  zed  Eigenvectors 

From  Eqs.(A&)  and  ( A24  )  it  follows  that 


CP  = 

p  J 

m2" 

UTP 

=  I 

i  A28 

pqt 

-  r 

■  42  Q 
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and 


T  -1 

Q  C  =  P  C 

=  P  ‘cpp  1 

=  J  0T  ■ A JO > 

Equation  (A30)  implies  that  Q  is  the  matrix  of  the  general  left 

eigenvectors,  while  Eqs.(A28)  and  (A29)  indicate  orthogonality  between  the 
right  and  left  general  i^ed  eigenvectors. 


Computation  of  General  1 zed  Right  and  Lef  t  E 1 qenvec  tor s 


Introducing  the  individual  generalized  eigenvectors  into  P  and  Q  suet 

that 

P  =  [  p1  p  2'.  .  .  p  L> .  .  .  .1  '  A 3 1  > 

and  Q=  t  q  l>  q  "’.  .  .  q‘  .  .  .  ]  (A32) 

it  can  be  demonstrated  that  the  application  of  Eqs.(A27)  and  (A30)  for  J  of 
the  Type  C"  t  i  i  J>  and  order  kxk  leads  to  the  following  relations: 


(C-AI  )p  =  X.p  or 

1  2 

(C-A.T  )p  =  Ap  or 

2  3 

1 C-X.I )  p  =  \p  or 

*k-i  *  k 

C-  I  )p  -  O  or 

T 

and  q  (OXl>  -  O  or 

i 

T  T 

q  (0X1  )  =  Xq  or 

*1  4 


(C-AI  )lp  -  O 
(C-AlX  *p  =  o 

7 

<  A33 ) 

(C-Xl)2p  =  o 

*k  1 

(C-AI  ip  •-  O 

*k 

q  X  C-AI  1  =  O 

qT i  C  -X  I  >2  =  O 
2 

( A3A  ) 


r  t 

q  ((’.■  /  I  )  -  q 
n  i  M  7 


q  'C-AI)  -■  A  q 
n  n  -i 


or  qT  <  C  ->  I  )  *  1  =  O 

*1  i 

or  qT<<:  I  '  =  O 


where  for  simplicity  all  super  sc r i d t s  t i >  have  been  omitted  and  the  subscripts 


?1 


1  through  k  refer  to  the  eigenvectors  comprising  the  columns  in  p  and  q 
Since  the  analysis  of  the  relative  strengths  in  combat  engagement s  invol.es 
the  calculation  of  the  dominant  left  eigenvector  q  it  is  clear  that  the 

i 

selection  of  J  as  the  lower  triangular  matrix  was  mo:  e  convex  1  ei  t  because  i  >. 
the  case  of  defective  matrice'  C  it  resulted  in  a  simple  e>p' ession 

q^C-Al  :  s  O 

or  (CT-.\  I  )q  =  O  , 

for  the  calculation  of  q  . 


d.  Ex  amp  1  es  of  the  Domi  rant  Le-f  t  E  1  gen-,  ec  tm  s 

The  dominant  left  e  1  qer.vec t or s  are  illustrated  here  for  se.eial  n  on.,, 
engagements  for  which  all  attrition  coefficients  are  equal  to  a  constant 
The  number  of  engagmq  t  vpes  of  weapons  is  indicated  b  .■  the  subscripts  ou.m 
with  each  matrix  C.  All  eigenvalues  of  T  are  shown,  with  the  dominant 
eigenvalue(s)  listed  first.  The  dominant  eigenvector  q  lias  been  norma  1 1 
on  its  largest  element  for  the  Blue  forces. 


C 


0  -b 

-  a  0 


=  /sb .  -/ati 


q  =  (  1  /h  a  > 

I  1 


0  0  -b 

-a  -a  (’ 


/Bab  •  -  /.Jaf  ; 


/Ph  a  ' 

3  , 1 


c 


0  0  .  .  .  -b 

0  0  . . . -b 


tJfi.l 


=  /mab 


•  '  i  /mat  ■ 


1  . 


-  Ymt  a 

•  l  l 


BE 


r  0 

0 

0 

0 

-b 

-b 

-b 

-b 

:  \  -  £  /  ab »  0*  0  ,  -2/a  b 

-  a 

-a 

0 

0 

} 

q  =  (  1  1  -  /b  /  a  -  /b  ya 

-a 

0 

0 

i 

1  ;  2 , 2 » 

1 

0 

0 

-b 

- 1 

c 

0 

0 

-b 

-b 

o 

0 

—  H 

-h 

"  a 

-a 

~  a 

0 

0 

i  ~a 

~  a 

-  a 

o 

0 

=  Vhab , 


q  =  ( 

1 


1  1 


0,  -  /<bab 

-  /3t77s7  - /3b, 


rs  i 


r.  - 

0 

0 

0 

-b 

-b 

-b 

;  ‘  1  /ab .  •"'«  0.  C,  -3/ab 

0 

c* 

o 

-  b 

-b 

-b 

o 

o 

0 

-b 

-b 

-b 

q  =  {  i  1  1  -  /b  /  a  -  /b  ,'a 

1  -a 

a 

-  a 

0 

o 

0 

1  a 

-  a 

-  a 

0 

0 

0 

! 

-a 

-a 

-  a 

o 

0 

0 

:  J.3> 

Tor  toe 

ger 

ier  a 

1  case 

of 

m  Blue  forces  or,  n  Red  forces 

coefficients  in  A  and  B 


>  J  , 


=  Vmnab 


1 1  -  (  1  i  .  im,.-,  /nt  na  /mb  na  .  ■  t.w.o-i 

i  -1 

The  following  two  examples  with  some  rero  attrition  coefficient1 
represent  cases  of  defective  matrices  C  with  repeated  eigenvalues. 


C 


0  0  -t>  -b  ;  V  -■  i/ab .  /at). 

0  0  0  -b 

a  0  0  C1  o  -  \  C>  1  0 

0  -a  0  0  s 


-  /at? .  -  /af> 
/  b  •  a  ; 

4x1 
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c 


0 

0 

0 

-b 

-b 

-b 

1  :  A  =  /ab,  /ah ,  0.  0, 

-  /at> 

0 

0 

0 

0 

0 

-b 

1 

1 

0 

0 

0 

0 

0 

-b 

q  =  {  0  0  1  0  0 

-  /b  ! 

•a 

0 

0 

0 

0 

0 

1 

0 

-a 

0 

0 

0 

0 

II 

0 

0 

-a 

0 

0 

0 

■;3,3/ 


e.  Compute  t  ion  of  the  Domnant  Let' t  Eigenvectors  and  tliejr  Eigenvalues 


T 

] n  addition  to  the  dominant  eigenvalue  >  the  matri  C  contains  also  an 

1 

eigenvalue  egual  to  - A  .  Consequently .  because  there  are  two  eigen. aloes  with 

equal  absolute  values  the  standard  power  method  of  computing  the  largest 

eigenvalue  and  its  corresponding  eigenvector  can  not  he  applied:  however,  the 

<,> 

power  method  with  a  shift  of  origin  car-  be  used  to  find  both  *•  and 

i  i 

Start  irq  with  Eq .  <C.35>  fo>*  i 

i 

(CT-X  1  1=0  '  L  .  io  1 

1 

and  adding  ul  to  Pott  sides,  where  u  is  a  constant  <e,q.«  1  >  , 


T 

L.  h. 


I  )  q  =  <  > 

i 


l:  q 


Hence , 


wher  e 


C1qj^qi 


t:T=o. 


\  -  -  u 

1 

In  this  way  the  shift  of  the  origin  for  thp  eigenvalues 
resulted  in  |  •.  +  u  |  being  greater  than  j  therefore.  the 

method  can  be  applied  to  Eq.  <C.38’  to  find  q  and  •  and  then 


>  E  .40  )  . 


and  hut 

a  i  v  1 1  i  q  ;  c  w  e 
t  •  (o. 
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FIGURES 

Fig.  1  Model  of  heterogeneous  force  engagements  ( Many-on-manv > . 

Fig.  2  Example  of  uncoupled  and  coupled  tactical  engagements. 

Fig.  3  Numerical  example:  one-on-two  engagement. 

Flu,  4  ‘iumei  ioul  example  (Cose  1 :  Variation  of  force  strengths  X  ,  <  .  X  , 

1  7.  3 

with  the  nondimensional  time  r  =  A.'^t  for  '£  =  l.fa  and  X  =  {  200  100  100  ). 

o  o 

Fig.  5  Numerical  example  fCase  2):  Variation  of  force  strengths  X  ,  X  ,  X  , 

with  the  nondimensional  time  r  =  A.'1  t  for  £  =  1.0  and  X  =  (  125  100  100  ). 

n  « 

Fig.  6  Numerical  example  (Case  3):  Variation  of  force  strengths  X  ,  x  ,  X  , 

12  3 

with  the  nondimensional  time  t  =  X'^t  for  ■£  =  0.8  and  X  =  {  100  100  100  '. 

o  o 

Fig.  7  Numerical  examples  (Case?.  1,  2,  and  3):  Variation  of  the  superior  its 
parameter  5  with  the  nondimensional  time  ;  =  X(1’t  for  -  1.6,  1.0,  and 

o.8. 
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a)  Uncoupled  engagements  Coupled  engagements 


Fig.  2  Examples  of  uncoupled  and  coupled  tactical  engagements. 
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Fig.  3  Numerical  example:  one-on-two  engagement. 
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Numerical  example  (Case  1):  variation  of  the  force  strengths  XA,  Xj,  and  X- 
vith  the  nondimensional  timeT  =>'^t  for  §)  =1.6  and  X  ={200  100  100}. 
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Numerical  example  (Case  2):  variation  of  the  force  strengths  a  ,  a  ,  a 
with  the  nondimensional  timer  =^'\  for  $0=l  O  and  X0={125  100  100}. 


Numerical  examples  (Cases  1,  2,  and  3):  variation  of  the  superiority 
parameter  S  with  the  nondimensional  timeT=$  for  $  =1.6,  1.0  and  0. 
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